AFOSR  ;TN-60-fi227 

■5  * 


AD433736 

CONTRACT  AP_ 49(638) -255 


( ARAP  Report  No.  51) 
EXAMINATION  OF  THE  SOLUTIONS 
OF  THE  NAVIER-STOKES  EQUATIONS 
FOR  A  CLASS  OF  THREE-DIMENSIONAL 
VORTICES  . 

Part  Ills  Temperature  Distribu¬ 
tions  for  Steady  Motion 


Coleman  duP.  Donaldson 
and 

Roger  D.  Sullivan  , 


Prepared  for 

Air  Force  Office  of  Scientific  Research 
Office  of  Aerospace  Research 
United  States  Air  Force 
Washington  25 >  D.  C. 

November  ”1963 

Aeronautical  Research  Associates  of  Princeton*  Inc. 
50  Washington  Road*  Princeton*  New  Jersey 


REPRODUCED  BY;  NT© 
U.S.  Department  of  Commerce 
National  Technical  Information  Service 
Springfield,  Virginia  22161 


ABSTRACT 


The  temperature  distributions  that  occur  inside  a 
steadily  rotated  cylindrical  vortex  tube  with  porous  walls 
having  steady  radial  inflow  of  an  essentially  incompressible 
fluid  are  obtained.  The  nature  of  the  distributions  for 
both  liquids  and  gases  is  discussed.  The  results  are  com¬ 
pared  with  Rott's  results  for  Burgers1  unbounded  vortex. 
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1 .  Introduction 

In  reference  1  (hereafter  referred  to  as  Part  l)  a 
class  of  solutions  of  the  Navier-Stokes  equations  for  the 
steady  flow  of  an  incompressible  viscous  medium  are  dis¬ 
cussed  in  some  detail.  These  flows  are  such  that  in  terms 
of  a  cylindrical  coordinate  system  ( r*  9S  z) *  the  respective 
velocity  components  (u*  v*  w)  are  of  the  form 

u  =  u(r) 

(1.1)  v  =  v(r) 

w  =  zw(r) 

and  the  pressure  is  of  the  form 

(1.2)  p  =  -|pCz2  +  p(r) 

where  p  is  the  density  and  C  is  a  constant.  It  has  been 
pointed  out  by  Donaldson  (reference  2)  that*  for  such  "essen¬ 
tially  incompressible"  flows*  the  energy  equation  can  be 
easily  solved  for  temperature  distributions  of  the  form 

(1.3)  T  =  T0(r)  +  z2T2(r) 

In  the  present  part  (Part  III)*  we  give  the  results  of 
a  brief  study  of  the  nature  of  such  temperature  distributions* 
partly  for  the  sake  of  completeness  and  partly  because  of  a 
general  interest  in  temperature  distributions  within  vortex 
flows.  The  results  that  are  obtained  are  closely  related  to 
the  work  of  Rott  (reference  3)  who  studied  the  temperature 
distributions  in  a  free  or  unbounded  vortex  of  the  form 
u  =  -ar*  v  =  v(r)*  and  w  =  2az.  The  present  study  shows 
that  the  temperature  distributions  in  confined  vortices  are 
very  similar  to  those  of  free  vortices. 


Attention  is  confined  to  vortices  of  the  family  desig 
nated  1-1  in  Part  X,  i.e.  single-cell  vortices  in  which 
the  flow  spirals  inward  from  a  porous  wall  at  r  =  R  and 
out  along  the  axis.  One  boundary  condition  on  the  tempera 
ture  is  taken  to  be 

(1.4)  T(R3  z)  =  T  =  constant 

The  only  other  condition  needed  is  that  the  temperature 
remain  finite  on  the  axis. 
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2.  The  Energy  Equation 

In  Part  I  (pp ,  16-22)  the  energy  equation  is  discussed 
in  detail.  Two  forms,  one  for  gases  and  one  for  liquids, 
which  are  consistent  with  the  constant  density  solutions  of 
the  continuity  and  momentum  equations  given  in  Part  I  are 
deduced.  The  equation  for  gases,  suitable  for  cases  in 
which  the  Mach  number  associated  with  the  tangential  veloc¬ 
ity,  v,  is  appreciable  while  that  associated  with  u  or  w 

*  .  . 

is  negligible,  is,  subject  to  (1.1), 


where  the  dissipation  is 


(2.2) 


As  usual  0^  denotes  the  specific  heat  at  constant  pressure 
and  k  the  thermal  conductivity,  both  assumed  constant. 

For  liquids  it  is  shown  In  Part  I  that  the  energy 
equation,  again  subject  to  (l.l),  is 


(2.3)  Po/u|^+w|^)  =  l/^|+i|E+i-f)+  •' 

P\  Sr  Sz /  \§r2  Sr  / 


*In  Part  I, 


should  read 


(1.5.23)  is  written  Incorrectly. 

+  +  Puv2 

St  r 

p 

S]D  _  puv1- 
St 


The  terms 


r 
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where  the  dissipation  is 


(2.4) 


As  pointed  out  by  Goldstein  (reference  4),  there  has 
been  some  confusion  in  the  literature  about  the  correct  form 
of  the  energy  equation  for  the  case  of  "nearly  incompressi¬ 
ble"  fluids.  The  forms  given  above  are  believed  to  be  the 
most  useful  and  are  consistent  with  both  the  recent  work  of 
Rott  (reference  3)  and  the  suggestions  of  Goldstein. 

The  necessity  of  assuming  that  the  Mach  number  associa¬ 
ted  with  u  or  w  is  negligible  leads  to  a  very  simple 
expression  for  the  dissipation  for  gases*  ©"*  which  makes 
the  solution  relatively  easy.  The  more  complicated  form 
for  liquids*  <pr*  makes  the  solution  more  involved  but 
gives  a  greater  richness  of  information. 

In  (2.1)  through  (2.4)*  u*  v*  and  w  are  known 
functions  that  are  determined  by  the  methods  outlined  in 
Part  I.  It  is  desirable  at  this  point  to  review  briefly  the 
way  the  solutions  of  family  1-1  are  obtained.  First*  the 
following  non-dimensional  variables  are  defined: 


(2.5) 

X  =  cr2/R2 

(2.6) 

F(X)  =  -  -|ur/v 

(2.7) 

H(X)  =  (vr/VR)H(c) 

where  V  is  the  value  of  v  at  the  wall*  v  =  q/p*  and 

(2.8)  c2  =  -CflV(l6v2) 
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Here  C  is  the  pressure-gradient  constant  appearing  in 
(1.2)5  it  is  always  negative  for  family  I-l*  so  c  is 
real.  The  equation  of  continuity  now  shows  that 

(2.9)  w  =  z w  =  4c(v/r)(z/R)f’ 

where  the  prime  denotes  differentiation  with  respect  to  X. 

As  shown  in  Part  I3  pp.  34-33*  F  and  H  are  governed 
by  the  equations 


(2.10) 

(XFU )  ’  +  FF”  -  F’2  +  1  =  0 

(2.11) 

+ 

ffl 

11 

0 

These  equations  are  integrated  numerically  with  the  initial 
conditions  given  by 

(2.12)  F(0)  =  0,  F’(0)  =  k3  F"(0)  =  A2  -  1 

(2.13)  H(0)  =  0S  H'(0)  =  B 

where  the  constant  A  is  (for  family  I-l)  between  0  and 

1,  and  the  constant  B  is  arbitrary.  Since  it  is  speci¬ 
fied  as  a  boundary  condition  that  w  =  0  at  the  wall  r  =  R, 
it  is  seen  from  (2.5)  and  (2.9)  that  F!(c)  =  0.  Hence  for 
each  A  in  (2.12)  integration  of  (2.10)  continues  to  the 
first  (for  family  I-l)  zero  of  Ff.  The  value  of  X  at 
which  it  occurs  determines  the  value  of  c  corresponding  to 
the  particular  A.  Integration  of  (2.1l)  can  be  performed 
simultaneously  or  separately.  The  other  characteristics  of 
the  flow  are  determined  from  the  knowledge  of  the  functions 
F  and  H. 

Accordingly  (2.1)  and  (2.3)  are  reivritten  in  terms  of 
the  dimensionless  variables  and  at  the  same  time  the  assump¬ 
tion  (1.3)  is  introduced  in  the  modified  form: 
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2 

(2.14)  T  =  T0(X)  +  Tp(X) 

U  R2  ^ 

It  is  found  that  each  term  is  either  independent  of  z  or 

p 

else  multiplied  by  (z/R)  .  Since  the  equality  must  hold  for 
all  values  of  z,  two  equations  are  deduced  from  each  of 
(2.1)  and  (2.3).  For  gases  they  are 


(2.15) 

(2.16) 


!  \  » 


(XTp 


Pr(FTg  -  2Fl12)  =  0 


(XTq)  '  +  PrFTo  -  -  i-  T£ 


+  Pr 


V2c 


opH2(o)L2  X2 


while  for  liquids  they  are 


2  2 

(2.17)  (XTg)*  +  Pr(FTp  -  2F!T2)  =  -l6  Pr  XF"2 

c  R 


(XT«)’  +  PrFTn  =  -  ~  Tn  -l6  Pr  -^-4 


L0 

(2.18) 


2c  2 


•p*2  -  f'JL  .  ( JL 

F  F  2X  +  V2X 


-  Pr 


V2c 


cpH2(c) 


(»'  -  ij' 


In  these  equations  Pr  is  the  Prandtl  number;,  |i.cp/k. 
In  both  sets  the  boundary  conditions  are  that  Tq 
and  Tg  be  finite  at  the  origin^  i.e. 


(2.19)  T0(0)  <  co,  T2(0)  <  co 

and.,  by  (1.4)  and  (2.5)s 

(2.20)  T2(c)  =  0;  T0(c)  =  Tw 
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3 .  Solution  in  the  Case  of  Gases 

We  define  the  function  G(x)  to  be  the  solution 
(determined  numerically)  of  the  differential  equation 

(3.1)  (XGr  )*  +  PrFG'  =  Prjf  ^  -  (V  -  £) 

with  the  initial  conditions 

(3.2)  G(0)  =  0,  G‘(0)  =  0 

Further,  we  define 


(3.3) 


A„  = 


«  (c ) 


G(x)  -  G(c) 


Then  it  can  be  verified  that  the  set 


(3.4)  T2  =  0 


(3.5) 


is  a  solution  of  (2.15)  and  (2.16)  which  satisfies  the  bound¬ 
ary  conditions  (2.19)  and  (2.20). 

The  computations  have  been  made  for  a  range  of  the 
Reynolds  number  =  UR/v,  where  U  =  u(R)  <0  (for  family 
I-l),  and  for  two  values  of  the  Prandtl  number,  Pr  =  l/2 
and  Pr  =  1.  Typical  distributions  of 


-  V 


v‘ 


(3.6) 


are  shown  in  Figures  3.1,  3.2 ,  and  3.3.  Figure  3.^-  shows 
the  variation  of  the  minimum  value 
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Ar 

min 


=  V  °> 


with  Reynolds  number.  These  results  are  discussed  in 
Section  5. 


Temperature  distribution  for  gases  at  =  -16.O. 


Fig.  3.3 
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4.  Solution  in  the  Case  of  Liquids 

The  technique  of  solution  in  the  case  of  liquids  is 
similar  to  that  just  outlined  for  gases,  but  considerably 
more  involved.  We  define  five  functions,  Ag,  Bg,  A^, 
AT,  and  Ajj,  through  their  differential  equations  and 
initial  conditions  as  follows 


(4.1)  (XAg) '  +  Pr(PAg  -  2F’Ag) 

(4.2)  Ag( 0)  =  C,  Ag(0) 

(4.3)  (XBg) '  -J-  Pr(FBg  -  2E,Bg) 

(4.4)  Bg( 0)  =  D,  Bg(0) 

(4.5)  (XAp*  +  PrPAp  =  -  Prj^F’2 

(4.6)  af(o)  =  0,  Ap( 0) 

(4.7)  (XA^)'  +  PrFA^  =  Bg 

(4.8)  At(0)  =  0,  A^(0) 

(4.9)  (XAg)’  +  PrFAy=  -  Pr  (e' 

(4.10)  Ah(0)  =  0,  A^(0) 


=  -PrXF"2 
=  2PrAC 

=  0 

=  2PrAD 


=  0 

like  B  in  (2.13),  may  be 


In  these  equations  C  and  D. 
given  any  value  which  is  convenient  in  the  numerical  Integra 
tion.  Further,  we  define 
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(4.11) 

A0  ~ 

16  c  . 

N2  1 
U 

af 

Ap(c)  + 

(4.12) 

A2  = 

16  c2 

Nu 

(a2 

-  “a) 

and 

(4.13) 

AV  = 

C 

l/H 

-  Ah(c) 

H2(c) 

where 

(4.14) 

K  = 

A2(c) 

Bp(c) 

Then  it  can  be  verified  that  the  set 


(4.15) 

(4.16) 


=  ^  +  ao  +  #av 

p  p 


is  a  solution  of  (2.17)  and  (2.18)  which  satisfies  the  bound¬ 
ary  conditions  (2.19)  and  (2.20). 

Again  the  computations  have  been  carried  out  for  a  range 
of  the  Reynolds  number  N^.  and  for  two  values  of  the  Prandtl 
number,  Pr  =  l/2  and  Pr  =  1.  Typical  distributions  of 
Aq,  A^,  and  A^.  are  shown  in  Figures  4.1,  4.2,  and  4.3. 
Figure  4.4  shows  the  variation  of  the  maxima  of  the  same 
quantities  with  Reynolds  number. 

From  (2.14),  (4.15),  and  (4.16)  we  find 


op(T 


-  T  )  = 
w' 


U‘ 


Ao  + 


(4.17) 


/ODJl 


Fig.  4.4  Maxima  of  the  A  functions  for  liquids  as  a 
function  of  Reynolds  number. 
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Thus  for  small  z/R,  Ay  may  be  Interpreted  as 

c  (T  -  T  ) 

(4.18)  Av=^L_wi 


If  V  is  sufficiently  large  compared  to  U  while  AQ  may 
be  interpreted  as 


(4.19) 


°P(T  -  V 

0  "  U2 


if  U  is  sufficiently  large  compared  to  V.  But  note  from 
Figures  4.1  and  4.4  that  when  J |  is  small  V  must  be 
very  much  larger  indeed  than  U  for  the  Ay  term  to  dominate, 
and  that  when  {N^j  is  large,  the  opposite  is  true.  Of 
course,  for  z/R  sufficiently  large,  Ag  may  be  interpreted 
as 


(4.20) 


.  M2  °ntT  -  V 

2  ~  \zj  n2 


These  results  are  discussed  further  in  the  following 
section. 
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5.  Discussion 

From  the  physical  point  of  view,  the  most  interesting 
aspect  of  the  results  is  probably  the  disparity  in  the 
behavior  of  liquids  and  gases.  That  is,  liquids  are  heated 
as  they  move  toward  the  center  but  gases  are  cooled.  This 
phenomenon  was  pointed  out  by  Rott  (reference  3)  for  a  free 
or  unbounded  vortex.  The  heating  in  the  case  of  liquids  is 
due  to  the  dissipation.  The  dissipation  is  still  present  in 
the  case  of  gases  of  course  but  it  is  more  than  counter¬ 
balanced  by  the  cooling  due  to  expansion  as  the  gases  approach 
the  low  pressure  region  near  the  center  of  the  vortex. 

In  reference  3  only  the  dissipation  due  to  the  tangen¬ 
tial  velocity  is  considered  in  the  solution  so  in  the  case  of 
liquids  it  is  appropriate  to  use  the  approximation  (4.l8)  in 
comparing  results.  But  A^.  in  that  approximation  and  A& 
for  gases,  by  (3.6)  represent  the  enthalpy  difference  divided 
by  V  ,  and  V  has  no  meaning  for  Burgers*  solution  (on 
which  Rott*s  work  is  based)  since  that  solution  represents  an 
unbounded  vortex.  Rott  presents  his  results  in  terms  of  the 
enthalpy  difference  divided  by  H*  where  H*  is  defined  in 
terms  of  the  circulation  at  infinity.  This,  of  course,  has 

no  meaning  for  the  present  solutions.  However,  Rott  shows 

p 

that  H*  =  1.22  vm  where  vm  is  the  maximum  value  of  v, 
and  this  provides  a  basis  for  comparison.  Accordingly,  the 
quantity  cp(T0  -  Tw)/(l.22  v^)  where  TQ  represents  the 
temperature  at  the  origin,  has  been  calculated  from  the 
formula 


for  gases,  and 


-  V 

1 

<1 

ro 

1.22  v2 
m 

1.22  v^ 
m 

%<T0  -  V 

i 

<1 

ro 

1.22  v2 

1.22  v£ 
m 

V°> 


v°> 
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for  liquids,  where  the  approximation  (4.18)  has  been  assumed. 
Figure  5.1  shows  the  variation  of  this  quantity  with  N^.. 

Also  indicated  in  the  figure  are  the  values  obtained  by  Rott 
for  liquids  with  infinite  Prandtl  number  and  for  gases  with 
Pr  =  l/2,  Pr  =  1,  and  Pr  =  <»,  These  are  labeled  =  -<» 
since  the  vortex  considered  has  no  boundary  at  a  finite 
radius . 

The  present  results  seem  to  approach  asymptotic  values 
as  j My  1  increases  but  it  is  not  clear  that  for  gases  they 
approach  the  values  given  by  Rott.  If  not,  the  difference  is 
probably  due  to  the  fact  that  at  large  Reynolds  numbers  the 
solutions  of  family  1-1  are  more  like  the  so-called  cosine 
solution  than  like  Burgers’  solution.  In  particular,  the 
axial  velocity  distribution  is  more  like 


w(z,  r) 
w(z,  0) 


cos 


than  is  like 


w(z,  r)  _  1 
w(z,  0) 

(See  Part  I,  p.  50  and  p.  64),  Thus  it  is  quite  likely  that 
the  temperature  difference  approaches  something  other  than 
Rott’s  value  as  approaches  minus  infinity. 

The  distributions  of  Ag  shown  In  Figures  4.1,  4.2,  and 
4.3  illustrate  an  interesting  effect,  namely  that  at  large 
enough  z,  by  the  approximation  (4.20),  the  maximum  tempera¬ 
ture  is  not  on  the  centerline  r  =  0  but  rather  at  a  radius 
greater  than  R/2.  This  behavior  is  due  to  the  fact  that  at 
large  z  the  fluid  entering  through  the  porous  wall  with 
w  =  0  must  turn  rapidly  and  acquire  a  substantial  axial 
velocity  component.  This  involves  a  large  contribution  to 
the  dissipation  near  the  wall.  In  fact  the  very  presence  of 
the  Ag  contribution  to  (4.17)  can  be  traced  to  the  term 
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p  p 

li(dw/dr)  (which  is  proportional  to  tt)  in  the  dissipation 
given  by  (2.4).  The  fluid  near  the  centerline  is  cooler 
since  it  entered  at  a  smaller  z  and  thus  experienced  a 
milder  dw/3r  and  less  dissipation. 

These  observations  serve  as  a  reminder  that  the  solu¬ 
tions  for  gases  are  not  valid  for  large  z  where  the  Mach 
number  associated  with  w  becomes  appreciable. 
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